Soliton and black hole solutions of su(iV) Einstein- Yang-Mills theory in anti-de Sitter 
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We present new soliton and hairy black hole solutions of su(N) Einstein- Yang-Mills theory in 
asymptotically anti-de Sitter space. These solutions are described by 7V+1 independent parameters, 
and have TV — 1 gauge field degrees of freedom. We examine the space of solutions in detail for su(3) 
and su(4) solitons and black holes. If the magnitude of the cosmological constant is sufficiently 
large, we find solutions where all the gauge field functions have no zeros. These solutions are of 
particular interest because we anticipate that at least some of them will be linearly stable. 
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I. INTRODUCTION 



Soliton and hairy black hole solutions of Einstein- 
Yang-Mills (EYM) theory and its variants have been the 
subject of extensive research since the discovery of non- 
trivial, spherically symmetric solitons [l| and 'colored' 
black holes in su(2) EYM in asymptotically fiat space- 
time. These black holes are 'hairy' in the sense that they 
have no magnetic charge, and are therefore indistinguish- 
able at infinity from a standard Schwarzschild black hole. 
There are discrete families of solutions, indexed by the 
event horizon radius T) x (with T) x = for solitons) and n, 
the number of zeros of the single gauge field function u>, 
each pair (r^ , n) identifying a solution of the held equa- 
tions. A key feature of the solutions is that n > 0, so that 
the gauge held function must have at least one zero (or 
'node'). These solutions, while they violate the 'letter' 
of the no-hair conjecture, may be thought of as not con- 
tradicting its 'spirit', since they are found to be unstable 
to classical, linear, spherically symmetric perturbations 
Q. There is also a large literature concerning analytic 
studies of the asymptotically flat su(2) EYM held equa- 
tions [3, H| , proving the existence of the above numerical 
solutions and detailed properties of the phase space of so- 
lutions. Since these initial discoveries a plethora of new 
soliton and black hole solutions have been found (see @ 
for a review). The present work combines two natural 
extensions of these initial studies: the generalization to 
su(N) EYM, and the inclusion of a negative cosmological 
constant. We now briefly review each of these general- 
izations in turn. 

Firstly, in asymptotically flat space, both charged and 
neutral numerical solutions of the su(N) EYM held equa- 
tions have been found [7( . We consider in this paper only 
purely magnetic solutions, which, in the asymptotically 
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flat case, are described by N — 1 gauge field functions 
uij. As in asymptotically flat su(2) EYM, solutions exist 
at discrete points in the parameter space, and can be in- 
dexed by the radius of the event horizon (if there is one) 
and the number of nodes of the ujj (all ujj having at least 
one zero). Once again, there is a general result [|[ that 
all these solutions must be unstable. The su(iV) EYM 
field equations are considerably more complicated than 
those for su(2) and correspondingly less analytic work 
has been done. Local existence of solutions of the field 
equations near the origin (for solitons), black hole event 
horizon (if there is one) and at infinity has been estab- 
lished 0, [l(| • There is a heuristic proof (following ) o f 
the existence of black hole solutions for general N 1JJ , 
but more rigorous work exits only for the case of su(3) 

The second generalization of asymptotically flat su(2) 
EYM that we consider in this paper is the inclusion of a 
non-zero cosmological constant A. When the cosmolog- 
ical constant is positive, soliton [l3| and black hole [3] 
su(2) solutions have been found. These solutions pos- 
sess a cosmological horizon and approach de Sitter space 
at infinity (for a complete classification of the possible 
space-time structures, see [15|). The phase space of solu- 
tions is again discrete, and the single gauge field function 
w must have at least one zero. Unsurprisingly, these solu- 
tions again turn out to be unstable [14l . [16[ . The inclusion 
of a negative cosmological constant (so that the space- 
time is asymptotically anti-de Sitter (adS)) may be moti- 
vated by recent progress in string theory, particularly the 
adS/CFT correspondence [I3|- It is found [3 EH that 
the solutions of su(2) EYM in adS possess quite differ- 
ent properties compared with their asymptotically flat or 
asymptotically de Sitter cousins. In particular, solutions 
for which the gauge field function uj has no zeros exist for 
sufficiently large |A|. Solutions exist in continuous open 
subsets of the parameter space, rather than at discrete 
points. In addition, for sufficiently large |A|, at least 
some of these solutions are stable under linear, spheri- 
cally symmetric perturbations [H, (this was subse- 
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quently extended to cover non-spherically symmetric lin- 
ear perturbations in [20j]). Therefore, while black holes 
cannot be given stable YM hair in either asymptotically 
flat or asymptotically de Sitter space, in asymptotically 
anti-de Sitter space, stable gauge field hair is possible. 

We are thus led to the following natural question: are 
there stable solutions of the S\x(N) EYM solutions with 
a negative cosmological constant? In this paper we will 
present the first soliton and hairy black hole solutions of 
su(N) EYM in adS, for N > 2. We consider only purely 
magnetic gauge fields, so that the YM field is described 
by N — 1 functions LVj. We make a detailed study of 
the phase space of solutions and their general properties 
in the particular cases N = 3, 4. Of particular interest 
is the existence, for sufficiently large |A|, of solutions in 
which all the ojj have no zeros. We anticipate that at 
least some of these solutions will be stable under linear, 
spherically symmetric perturbations, and will examine 
their stability in detail elsewhere. 

The outline of this paper is as follows. In section 
HIl we discuss the field equations, our ansatze for the 
fields and the boundary conditions that must be sat- 
isfied, considering the cases of black holes and solitons 
separately. Our new solutions are discussed in detail in 
section IIII1 focussing particularly on the phase space of 
solutions. Finally, our conclusions can be found in sec- 
tion llVI Throughout this paper, the metric has signature 
(— , +, +, +) and we use units in which AnG = c = 1. 



II. ANSATZ, FIELD EQUATIONS AND 
BOUNDARY CONDITIONS 



A. Ansatz and field equations 

We consider static, spherically symmetric, four- 
dimensional solitons and black holes with metric 



ds 2 



liS A dt 2 + fi- 1 dr 2 + r 2 de 2 + r l sin 2 9 d<j) z , (2.1) 



where the metric functions fi and S depend on the radial 
co-ordinate r only. In the presence of a negative cosmo- 
logical constant A < 0, we write the metric function /i 
as 



2m(r) Ar 2 



(2.2) 



The most general, spherically symmetric, ansatz for the 
su(N) gauge potential is plj ]: 

A = Adt + Bdr+^(C-C H ) dO 

- % - [(C + C H )sm9 + Dcos8] d<p, (2.3) 

where A, B, C and D are all (N x N) matrices and C H 
is the Hermitian conjugate of C . The matrices A and 
B are purely imaginary, diagonal, traceless and depend 



only on the radial co-ordinate r. The matrix C is upper- 
triangular, with non-zero entries only immediately above 
the diagonal: 

C7 ili+1 = Wj -(r) e ^W, (2.4) 

for j = 1, . . . , N — 1. In addition, I? is a constant matrix: 

D = Diag (N - 1, N - 3, . . . , -N + 3, -N + 1) . (2.5) 

Here we consider only purely magnetic solutions, so we 
set .4 = 0. We may also take B = by a choice of 
gauge [2l| . From now on we will assume that all the 
LOj{r) are non-zero (see, for example, Q for the possibil- 
ities in asymptotically flat space if this assumption does 
not hold). In this case one of the Yang-Mills equations 
becomes UM 







Vj 



1. 



,JV-1. 



(2.6) 



Our ansatz for the Yang- Mills potential therefore reduces 
to 



A = ^(C -C H ) d6- l -[(C + C H ) sin(9 + £>cos( 



(2.7) 



where the only non-zero entries of the matrix C are 

C jtj+1 = Wi (r). (2.8) 

The gauge field is therefore described by the N — 1 func- 
tions Wj(r). We comment that our ansatz (|2.7[) is by 
no means the only possible choice in su(N) EYM. Tech- 
niques for finding all spherically symmetric su(N) gauge 
potentials can be found in [22| , where all irreducible mod- 
els are explicitly listed for N < 6. 

With the ansatz (|2.7p . there are N— 1 non-trivial Yang- 
Mills equations for the N — 1 functions u>j : 



r 2 ^ 



2m - 2r 3 pg 



2Ar 3 \ 

—ju'. + Ww^O (2.9) 



for j = 1, . . . , N — 1, where a prime ' denotes d/dr, 



Po 



1 N 
47^ [ 



' 2 2 



N-1 + 2j) ,(2.10) 



i=i 



Wj - 1 



1 



and ojq = ljn = 0- The Einstein equations take the form 

m = fiG + r p 8 , — = — , 
b r 

where 

iV-l 



(2.12) 



rl 



(2.13) 



j'=i 



Altogether, then, we have N + 1 ordinary differential 
equations for the N + 1 unknown functions m[r), S(r) 
and u)j(r). 
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The field equations (]2.9I2.12[) are invariant under the 
transformation 

Wi (r) - - Wi (r) (2.14) 
for each j independently, and also under the substitution: 
3->N-j. (2.15) 

B. Boundary conditions 



The field equations (|2. 912.12ft are singular at the origin 
r = (for regular, soliton solutions), the black hole event 
horizon r = Th (if there is one) and at infinity r — > oo. 
We therefore now discuss the boundary conditions that 
must be satisfied by the field variables at these singular 
points. Local existence of solutions of the field equations 
in neighborhoods of these singular points will be rigor- 
ously proved elsewhere [23[, generalizing the local exis- 
tence proofs in the asymptotically flat case 0, |T(| ■ The 
boundary conditions satisfied by black hole solutions are 
more easily stated, so we consider those first. 



1. Black holes 

We assume there is a regular, non-extremal, black hole 
event horizon at r = rh, where fi(r) has a single zero. 
This fixes the value of m(r^) to be: 



2m(r h ) = r h ^. 



(2.16) 



The field variables uj j(r), m(r) and S(r) will have regular 
Taylor series expansions about r = r^: 

m(r) = m(r h ) + to' \r h ) (r ~ r h ) + O (r - r h ) 2 ; 
Uj{r) = ^j{r h ) + 0{r-r h ) ; 

S{r) = S{r h ) + 0{r-r h ). (2.17) 

Setting jtx(rh) = in the Yang-Mills equations (|2.9[) fixes 
the derivatives of the gauge field functions at the horizon: 



Wjir^ujjirh) 



2m{r h ) - 2r\p e {r h ) 



2Ar? 



(2.18) 



Therefore the expansions (|2.17|) are determined by the 
N+l quantities ujj(rh), rh, S(rh) for fixed cosmological 
constant A. For the event horizon to be non-extremal, it 
must be the case that 



2m'(r h )=2r 2 h pg(r h ) < 1 - Ar 2 h , 



(2.19) 



which weakly constrains the possible values of the gauge 
field functions Uj(rh) at the event horizon. Since the field 
equations (|2.9I2.12| are invariant under the transforma- 
tion (|2.14|) , we may consider wArh) > without loss of 
generality. 



At infinity, the boundary conditions are considerably 
less stringent than in the asymptotically flat case. In 
order for the metric (|2.ip to be asymptotically adS, we 
simply require that the field variables Wj(r), m(r) and 
S(r) converge to constant values as r — > oo, and have 
regular Taylor series expansions in r~ 1 near infinity: 

m(r) = M + O^- 1 ); S(r) = 1 + O (r^ 1 ) ; 
Uj(r) = wj.oo + O (r- 1 ) . (2.20) 

Since A < 0, there is no cosmological horizon. 



2. Solitons 

Soliton solutions have the same boundary conditions 
(|2.20p as r — > oo as black hole solutions. The boundary 
conditions at a regular origin, however, are considerably 
more complicated than at a black hole event horizon or 
at infinity. For the asymptotically flat case, they have 
been derived in @. As may be expected, the modifica- 
tions required by the presence of a non-zero cosmological 
constant are not great. However, given the complexity of 
these boundary conditions, we now describe their deriva- 
tion in some detail. 

We begin by assuming a regular Taylor series expan- 
sion for all field variables near r = 0: 

m(r) = rriQ + mir + rri2r 2 + 0(r 3 ); 
S(r) = S + Sir + S 2 r 2 + 0(r 3 ); 
Wj(r) = Ujfl + u jtl r + Uj, 2 r 2 + 0(r 3 ); (2.21) 

where the rrn, Si and tUj^ are constants. The expansions 
(|2.21[) are substituted into the field equations (|2. 912. 12|) 
to determine the values of the constants. The constant 
So is non-zero in order for the metric to be regular at the 
origin, but otherwise arbitrary since the field equations 
involve only derivatives of S. 

Regularity of the metric and curvature at the origin 
immediately gives: 

toq = mi = ?7i2 = 0, Si — 0, ujj t i = (2.22) 



and 



,-,0 = ±y/j(N-j). 



(2.23) 



Without loss of generality (due to Q2.14[l ). we take the 
positive square root in (|2.23[) . 

Examination of the leading order terms in the Yang- 
Mills equations (|2.9p gives the following constraint on 

UJ2 = (^1,2, ■ ■ ■ ,^n-i,2) T - 

M N -ilj 2 = 2uj 2 . (2.24) 
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Here, Mn-i is the (N — 1) x (N — 1) matrix 



M 



N-l 



2(N-1) -y/(N- 1)2(N- 2) • 

y/(N — 1)2(N- 2) 2.2(iV - 2) -^2(N - 2)3(Y - 3) • 

--^2(^-2)3(^-3) 2.3(JV-3) 



V 



o 





■y/(N-l)2(N- 

2(N-1) 



Therefore u> 2 is an eigenvector of the matrix Mn-i with eigenvalue 2 if one exists, otherwise u) 2 = 0. 
To find the eigenvalues and eigenvectors of .Mjv-ij we nrs t note that it can be written in the form 



where 



V 



M 



N-l 



N-l 



M N -i = V N -xM N -xDx\, 



Diag (VN~\, y/2(N-2), y/3(N-3), 

( 2(N - 1) -2(JV-2) 
-(N-l) 2.2(JV-2) -3(JV-3) 
-2(JV-2) 2.3(JV-3) 



-(N-l) 
2(N-1) J 



(2.25) 



(2.26) 



(2.27) 



Then the matrices Mn-i and -Mn-i have the same 
eigenvalues, and the eigenvectors of Mn-i can be de- 
duced from those of M. jv-i- This result is useful because 
the matrix M. n-\ has been studied in detail in [§]. There 
it is proved that the eigenvalues of M n-i (and therefore 
those of Mn-i) are: 

8 k = k(k + l), k = l,...,N-l. (2.28) 

The eigenvectors of Mn-i in general involve Hahn poly- 
nomials [24j . and can be found explicitly in [9J. The 
eigenvectors v k for N — 3, 4 will be presented in sections 
IIII C 11 and IIII C 21 when we discuss the soliton solutions 
for su(3) and su(4) EYM, respectively. 
Therefore, we set 

w 2 = Mi, (2-29) 

where V\ is a (suitably normalized) eigenvector of A^jv_i 
with eigenvalue 2, and b\ is an arbitrary constant. From 
the Einstein equations (|2.12p , we find that 7713 and £2 are 
fixed and given in terms of the ojj^- 

Expanding the gauge field functions u>j to order r 2 
has therefore only introduced one arbitrary parameter, 
namely b\. However, it is expected that N — 1 indepen- 
dent parameters will be required to describe the N — 1 
independent functions ojj(r). Therefore, we must work 
to higher order in r in order to introduce more arbitrary 
parameters. 

Considering the next order in the Yang-Mills equations 
(|2.9|) . and setting 0)3 = (^1,3, ■ ■ • ,ojn-i,3) T , we find 



so that we may set 

u 3 = b 2 v 2 , (2.31) 

where 62 is an arbitrary constant and v 2 is an eigenvector 
of Mn-i with eigenvalue 6 (k = 2 in ([2~28|D . The Ein- 
stein equations (|2.12p are then used to determine 1714 and 
S3 (which will also depend on the cosmological constant 
A)._ 

Since we require N — 1 arbitrary parameters for the 
N — 1 independent functions Wj(r), the above analysis 
therefore suggests that we need to expand the ojj(r) up 
to r N+1 in order to have N — 1 arbitrary parameters in 
the expansion. This turns out to be the case, and a de- 
tailed proof will be given elsewhere [23[ . Determining the 

ujk = (u\,ki ■ ■ ■ , UN-i,k) T for k > 3 is slightly more com- 
plicated than for fc = 2,3 as outlined above. Examining 
the Yang-Mills equation (|2.9|) to fcth order, we find an 
equation for the u)k+i of the following form 

[M N -i-k(k + l)]u k+1 =c k+l (2.32) 

where c^+i is a complicated vector depending on 
o>i, . . . , u>k and 7713, . . . , m^. For A = 0, the form of Ck+i 
is given explicitly in [9]; when A < there are minor 
modifications which we do not write here (they will be 
given in [HI]). Since, in later sections, we present solu- 
tions just for su(3) and su(4) EYM, we will not need the 
exact form of the cjt+i. If Vk is an eigenvector of VWjv_i 
with eigenvalue £k (|2.28[) . we can solve equation (|2.32p 
for Wfc+i: 



Mn-i^3 = 6o; 3 , 



(2.30) 



Wfc+i = b k v k + Uk+i, 



(2.33) 
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where Uk+i is a particular solution of (|2.32[) chosen by re- 
quiring that Mfc+i is a linear combination of V\ 1 . . . , ffc-i- 
It is proven in Q that there is a unique solution of (|2.32|) 
subject to this constraint, in the A = case. This can 
be extended to A < 0, but we do not present the lengthy 
details here (23l |. 

The upshot of all this analysis is that the expansion 
of the fields, near the origin, is written as follows (where 
CO = (wi, . . . ,Wjv-i) ): 



m(r) 
S(r) 



m 3 r 3 + 0{r 4 ); 
So + O(r 2 ); 

N-l 



w(r) = u> + J2 hv k r k+1 + 0(r N+1 ), (2.34) 



fc=i 



where 



w 



(ViV^T, v/2(7V - 2), . . . , y/{N - 1)) T . (2.35) 



The expansions (|2.34|) give the field variables in terms of 
the N — 1 parameters b\, . . . , &/v— l an( i are those which 
are used in the numerical integration of the field equa- 
tions in the next section. 



III. SOLUTIONS 



The field equations (|2.9I2.12[) have the following trivial 
solutions. Setting cjj (r) = ±y/j(N — j) for all j gives the 
Schwarzschild-adS black hole with m(r) = M = constant 
(which can be set to zero to give pure adS space). Set- 
ting uij (r) = for all j gives the Reissner-Nordstrom-adS 
black hole with metric function 



2M Q At 2 



and magnetic charge Q given by 

Q 2 = -N(N+l){N-l). 
6 



(3.1) 



(3.2) 



There is an additional special class of solutions, given by 
setting 



(r)=±y/j(N-j)oj(r) Vj = l,...,N-l. (3.3) 



In this case, we follow [9( and define 



\n = ^N(N-1)(N + 1), 
and then rescale the field variables as follows: 



R = A^V; A = X 2 N A; 



(3.4) 



W 1 ' 

Cj(R) = u{r). 



m{R) = A^m(r); S(R) = S(r); 



Note that we rescale the cosmological constant A (this is 
not necessary in Q as there A = 0). The field equations 
satisfied by m(R), S(R) and u)(R) are then 



dm 
~dR 



uG + R 



. 2 d 2 UJ 



<2~ 



° = R ^dR 2 



1 dS _ 
SdR ~ 

2m - 2R i p 9 - 



2G 
~~R ; 
2Ai? 3 



dR 



+ \1-Cb 2 ] Q, 



where we now have 

H = l 



2m 



AR 2 



and 



G 



duj 
dR 



2R? 



- UJ 



(3.6) 



(3.7) 



(3.8) 



(3.5) 



The equations (|3.6p are precisely the su(2) EYM 
field equations. Furthermore, the boundary conditions 
(|2. 1712.2012. 34p also become those for the su(2) case. 
This is straightforward to see for the boundary conditions 
at the horizon (|2.17p or at infinity (|2.20|) . At the origin 
(|2.34p . the su(2) embedded solutions are given by b\ ^ 0, 
but 62 = • • • = &JV-1 = 0. Therefore any su(2), asymp- 
totically adS, EYM soliton or black hole solution can be 
embedded into su(7V) EYM to give another asymptoti- 
cally adS soliton or black hole. We will see later in this 
section how the embedded su(2) solutions fit in the solu- 
tion spaces for larger N. 

To find genuinely su(N) solutions, the field equa- 
tions (|2.9l2.12p are integrated numerically using standard 
'shooting' techniques [251 ] . The equation for S(r) decou- 
ples from the other Einstein equation and the Yang-Mills 
equations so can be integrated separately if required. For 
su(N) solutions, we therefore have N coupled ordinary 
differential equations to integrate (N—l Yang-Mills equa- 
tions and one Einstein equation). For black holes, we 
start integrating just outside the event horizon, using as 
our shooting parameters the N variables Wj(r^), r/i, sub- 
ject to the weak constraint (|2.19[) . For solitons, we start 
integrating close to the origin, using as our shooting pa- 
rameters the (N — 1) variables bj (|2.34[) . In the soliton 
case, there are no a priori bounds on the parameters 
bj. In cither case, the field equations are then integrated 
outwards in the radial co-ordinate r until either the field 
variables start to diverge or they have converged to the 
asymptotic form at infinity. 

We now turn to a detailed discussion of the solutions 
we find. As well as presenting some examples of solutions, 
our particular focus in this section will be the structure 
of the space of solutions, as a subset of the phase space of 
parameters characterizing the solutions. We will examine 
the solution spaces in detail for su(3) and su(4) solitons 
and black holes, focusing on the numbers of zeros of the 
gauge field functions. The solution spaces we present 
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may not necessarily be complete, as our approach has 
been to scan the parameter space using a grid. It is 
therefore possible that solutions in which the gauge field 
functions have different numbers of zeros exist between 
the points of our grid. However, our figures will reveal 
the key features of the solution spaces. Of particular 
interest will be the existence of solutions where all the 
gauge field functions have no zeros. 

A. su(2) solutions 

We begin by reviewing the phase space of su(2) solu- 
tions, in which case we have a single gauge field function 
uj(r). Many of the properties we find in the phase space of 
solutions for su(N), N > 2 are also seen in the su(2) case 
and so it is informative to examine this simpler situation 
first. 

1. su(2) solitons 

Near the origin, one parameter, b, is required, and the 
expansion (|2.34[) reduces to 

uj(r) = 1 + br 2 + O (r 3 ) . (3.9) 

For solitons, the phase space has been studied in detail by 
[26| . We have verified their results and the phase space 
is shown in figure [1] (note that our parameter b in (|3.9p is 
equal to —b in [26l]). The phase space is parameterized by 
just two quantities: the cosmological constant A and the 
shooting parameter b (13. 9[) . The shaded regions in figure 
[T] indicate those values of the parameters for which we 
were unable to find a regular solution all the way out to 
infinity. Where we do find solutions, they occur in open 
subsets of the plane. We label these open sets by n, the 
number of zeros of the single gauge function to. We draw 
the reader's attention to the following particular features 
of the soliton phase space: 

1. The number of zeros of the gauge field function 
increases as |A| decreases or b decreases. 

2. Solutions in which to has no zeros occur for suffi- 
ciently large |A|. 

3. As | A| decreases, we find fewer solutions. The phase 
space breaks up into smaller and smaller regions. In 
the limit A — > 0, we are left with solutions just at 
discrete points, which are the Bartnik-McKinnon 
solitons in asymptotically flat space [l|. 

Note that the solution with 6 = exists for all A, and 
simply corresponds to pure adS, with Lu(r) = 1. 

2. su(2) black holes 

We next turn to the phase space of su(2) black hole 
solutions. There are now three parameters describing 



-1 
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n = 2 




n = 
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log 10 {-A) 

FIG. 1: Phase space of soliton solutions of su(2) EYM. The 
shaded area denotes those pairs (A, b) (where b is the shooting 
parameter giving the form of the gauge field function u> near 
the origin) for which no regular solution is found. The un- 
shaded regions correspond to regular solutions, the number of 
nodes n of the gauge field function lu being indicated for each 
region. For values of b just below the region in which n = 3 
we found solutions for which n = 4, but the latter region is 
too small to be seen on the graph. 

the solutions, r^, A and w(r^). In order to plot two- 
dimensional figures, we fix either r/j or A and vary the 
other two quantities. For su(2) black holes, the constraint 
(|2.19p on the value of the gauge field function at the event 
horizon reads 

(^rtf-lf <rl(l-Arl). (3.10) 

Whether we are varying or A, we perform a scan over 
all values of ojh which satisfy (|3.10[) . 

Firstly, we show in figure [5] the space of black hole 
solutions for fixed A = —0.01 and varying event horizon 
radius r^. The outermost curves in figureOare where the 
inequality (|3 . 1 0|) is saturated. Immediately inside these 
curves we have a shaded region, which represents values 
of (rh, w(rh)) for which the constraint (|3.10[) is satisfied, 
but for which we are unable to find black hole solutions 
which remain regular all the way out to infinity. As with 
the solitons, where we do find solutions, we indicate in 
figure [5] the number of zeros of the gauge field function 
u>(r). The solution for which uj(rh) — 1 is simply the 
Schwarzschild-adS black hole, while that for u>(rh) = is 
the magnetically charged Reissner-Nordstrom-adS black 
hole, as described above. The following key features are 
apparent from figure O 

1. We find solutions in which the gauge field function 
has more zeros as we decrease 77, or u)(rh)- 

2. As rh — > 0, the constraint (|3.10|) implies that 
w(r/i) — > 1, as can be seen in figure [H This is 
because the black hole solutions become solitons in 
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./no solution 






n - 2 





n = 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

lo 9lo( r h) 




1 2 

iog 10 (-A) 



FIG. 2: The space of su(2) black hole solutions when A = 
—0.01, for varying ru- The shaded region indicates values of 
the gauge field function ui(rh) at the event horizon for which 
the constraint (|3.10|) is satisfied, but for which we find no 
well-behaved black hole solution. The number of zeros n of 
the gauge field function uj are indicated in those regions of the 
phase space where we find black hole solutions. Elsewhere on 
the diagram, the constraint (|3.10|l is not satisfied. Between 
the region where n — 2 and the shaded region we find black 
hole solutions with n — 3, 4 and 5, but these regions are too 
small to indicate on the graph. 



FIG. 3: Phase space of su(2) black holes with ru = 1 and vary- 
ing A. The shaded region indicates values of the gauge field 
function U)(rh) at the event horizon for which the constraint 
(|3.10p is satisfied, but for which we find no well-behaved black 
hole solution. The number of zeros n of the gauge field func- 
tion uj are indicated in those regions of the phase space where 
we find black hole solutions. Elsewhere on the diagram, the 
constraint (|3.10[) is not satisfied. As well as the regions where 
n = 0, . . . , 4 as marked on the diagram, we find a small region 
in the bottom left of the plot where n — 5. This region is too 
small to indicate on the current figure, but can be seen in 
figure g] 



this limit, and, for solitons, we have w(0) = 1 (|3.9p . 
However, as can be seen in figure [TJ for this value 
of A, there are different soliton solutions, with u> 
having different numbers of zeros. This feature is 
not readily apparent from figure [2] 

3. The phase space of solutions breaks up into smaller 
regions as Th decreases. 

We find similar behaviour on varying for different val- 
ues of A. 

We now fix the event horizon radius = 1 and vary 
the cosmological constant A. The solution space in this 
case is shown in figure^ with a close-up for smaller values 
of |A| in figure 2] Once again, in figures [3] and 0] we 
have shaded those regions where the constraint (|3.10p 
is satisfied, but no regular black hole solutions could be 
found. Where we do find solutions, the number of zeros 
of the gauge field function u(r) is indicated in the figures. 
Similar behaviour is observed as in the soliton case (figure 
[1]), namely: 

1. The number of zeros of the gauge field function 
increases as |A| or w(r^) decreases. 

2. As A — > 0, the phase space breaks up into discrete 
points, which correspond to the asymptotically flat 
'colored' su(2) black holes @. 




n = 



4-2 

iog 10 (-A) 



FIG. 4: Close-up of the phase space of su(2) black holes with 
Vh, = 1 and smaller values of A. In the bottom left of the plot 
there is a small region of solutions for which n = 7, but the 
region is too small to be visible. 



3. For sufficiently large |A|, we find solutions in which 
the gauge field function has no zeros. 
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FIG. 5: Typical su(3) black hole solution, with = 1, A = 

-1, u/i(r h ) = 1.2 and uj 2 {r h ) = 1.3. In this example, both FIG- 6: Example of an su(3) black hole solution, with r h = 1, 
gauge field functions have no zeros. A = -0.0001, u>i(r h ) = 1.184 and Lo 2 {r h ) = 1.216. In this 

case, both gauge field functions have three zeros. 



B. Black holes 

We now turn to solutions of the su(N) EYM field equa- 
tions with N > 2, considering firstly black holes and then 
solitons. Many of the features of the su(2) solutions out- 
lined in the previous section will be replicated for larger 
N. 



1. su(3) black holes 

For 5u(3) EYM, there are two gauge field functions 
LO\{r) and W2(V)j and therefore four parameters describ- 
ing black hole solutions: r^, A, wi(r/,) and ^(rh)- Us- 
ing the symmetry of the field equations (|2.14[) . we set 
LUi(rh),u)2{rh) > without loss of generality. The con- 
straint (|2.19[) on the values of the gauge field functions 
at the horizon becomes, in this case: 

[cu 1 (r h ) 2 -2] 2 +[uj 1 (r h ) 2 -uj 2 (r h ) 2 } 2 
+ [2-cu 2 (r h ) 2 ] 2 

<2r 2 (l~Ar 2 ). (3.11) 

Two typical black hole solutions are shown in figures [5] 
and[6l The metric functions behave in a very similar way 
to the su(2) solutions [lfl [l9j], smoothly interpolating 
between their values at the horizon and at infinity. We 
note that S(r) in particular converges very rapidly to 
1 as r — > oo. In figure [5l we show an example of a 
black hole solution in which both gauge field functions 
have no zeros. We note that both gauge field functions 
are monotonic, however, one is monotonically increasing 
and the other monotonically decreasing. In our second 
example (figure [5]) both gauge field functions have three 
zeros. Although, in both our examples the two gauge 
field functions have the same number of zeros, we also 



find solutions where the two gauge field functions have 
different numbers of zeros (see figures [8] and [9|) . 

We now examine the space of black hole solutions. 
Since we have four parameters, in order to produce two- 
dimensional figures, we need to fix two parameters in 
each case. We find that varying the event horizon radius 
produces similar behaviour to the su(2) case, so for the 
remainder of this section we fix rh — 1 and consider the 
phase space for different, fixed values of A, scanning all 
values of wj.(r/,), u>2(rh) such that the constraint (|3.11[) is 
satisfied. From the discussion at the beginning of section 
IIII1 we have embedded su(2) black hole solutions when, 
from (pO| : 

u) X (r) = V2u(r) =u 2 {r) (3.12) 

which occurs when cji(r/j) = uJ%{rh)- 

In figures [TlfTUl we plot the phase space of solutions 
for fixed event horizon radius r/, = 1 and varying cos- 
mological constant A = —0.0001, —0.1, —1 and —5 re- 
spectively. In each of figures [7]fT0l we plot the dashed 
line uo\(rh) = cj 2 (f/i); along which lie the embedded su(2) 
black holes. It is seen in all these figures that the solution 
space is symmetric about this line, as would be expected 
from the symmetry (|2.15[) of the field equations. 

As in the su(2) case, for small values of A (see figure 
[7J the solution space fragments and we find very few so- 
lutions. The values of (wi(r^), LU2{fh)) for which we find 
regular black hole solutions are indicated by black dots 
in figure [7] Above the main group of solutions, there 
can clearly be seen a couple of smaller regions of solu- 
tions. There is also a small region centered on and very 
close to the dashed line at about wi(r^) ~ 1.27, and the 
Schwarzschild-adS solution at cJi(rft) = Wa^/i) = V2- 
For this value of A, we find very complicated behaviour 
in the numbers of zeros (nj., n%) of the gauge field func- 
tions wi(r), u>2(r) respectively. We have found at least 
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FIG. 7: Solution space for su(3) black holes with = 1 
and A = —0.0001. The dashed line indicates where 0Ji{r%) = 
wj(r^), along which lie the embedded su(2) solutions. The 
black regions indicate where we have regular black hole solu- 
tions; elsewhere we find no solutions. In this case we find a 
wide variety of numbers of zeros of the gauge field functions, 
and so do not indicate all the different possibilities. For these 
values of A and th we find no solutions for which the gauge 
field functions have no zeros. The key feature in this figure 
is the fragmentation of the solution space and the fact that 
there are comparatively few solutions. 



FIG. 8: Solution space for su(3) black holes with r;, = 1 and 
A = —0.1. The numbers of zeros of the gauge field functions 
for the various regions of the solution space are shown. For 
other values of u>%(rh), wzirh) we find no solutions. There is 
a very small region containing solutions in which both gauge 
field functions have no zeros, in the top-right-hand corner of 
the plot. 



fourteen different combinations of the numbers of zeros 
of the gauge field functions, some of which occur only 
in very small regions of the parameter space. This be- 
haviour is too complicated to depict accurately in figure 
[7J The numbers of zeros of the gauge field functions vary 
between 1 and 4 (we find no solutions in which either 
gauge field function has no zeros). We stress that the 
gauge field functions do not have to have the same num- 
bers of zeros, for A = —0.0001 wc find that \n\ — r^l 
varies between and 2. 

The solution space is found to be symmetric about the 
line uji(rh) = ^(fh) n ot only in terms of where we find 
solutions, but also in terms of the numbers of zeros of the 
gauge field functions. To state this precisely, suppose 
that at the point <jj\(rh) = ffli, ^(fh) = a 2 we find a 
black hole solution in which u>\(r) has n\ zeros and u>2(r) 
has Tii zeros. Then, at the point cJi(r) = 02, u>2(f) = a ii 
we find a black hole solution in which uj\ (r) has zeros 
and Ui(r) has ri\ zeros. This is clearly seen in figures [8] 
and[5J and follows from the symmetry (|2.15p of the field 
equations. 

As we increase |A|, we find (see figures l8l fl"0|) that the 
solution space expands as a proportion of the space of 
values of uii(rh), ^(r/J satisfying the constraint (|3.11[) . 
It can also be seen from figures I51TTU1 that the number 
of nodes of the gauge field functions decreases as |A| in- 
creases, and that the space of solutions becomes simpler. 

For A = —0.1, there is a very small region of the so- 
lution space where both gauge field functions have no 





no solution 






n 1 = 1, 








n 2 = 

n i = ^<^<t 


*\T / / / 
^ / y .y 

/ V / 


/ / 


' n, = 0, 
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( r h) = «2( r h) 


\n 1 = 0, 
n 2 = 1 





0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



FIG. 9: Solution space for su(3) black holes with ru = 1 and 
A = — 1. The shaded region indicates where the constraint 
(|3.1ip is satisfied but we do not find black hole solutions. 
Outside the shaded region the constraint (|3.11|l does not hold. 
Where there are solutions, we have indicated the numbers of 
zeros of the gauge field functions within the different regions. 
For this value of A there is a large region in which both gauge 
field functions have no zeros. 



zeros. This region expands as we increase |A|, until for 
A = — 5, both gauge field functions have no zeros for all 
the solutions we find. 
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FIG. 10: Solution space for su(3) black holes with ru = 1 
and A = —5. It can be seen that for the vast majority of the 
phase space for which the constraint (|3.11[) is satisfied, we 
have black hole solutions in which both gauge field functions 
have no zeros. 



2. su(4) black holes 



FIG. 11: Solution space for su(4) EYM black holes with r h = 
1, A = —1 and wi(rfe) = 1.6. Where there are solutions, 
the numbers of zeros of the three gauge field functions are 
indicated for the relevant regions. Elsewhere in the figure we 
find no black hole solutions. As well as the regions indicated, 
we also find small regions where the numbers of zeros of the 
gauge field functions are (711,712,713) = (1, 1,0) and (0, 1, 1). 



In this case there are three gauge field functions and so 
the parameter space is five-dimensional. The constraint 
(|2.19[) satisfied at the horizon by the gauge field functions 
now reads: 

[ Wl (r,) 2 - 3] 2 + [u 2 {r h f ~ uj^ruf - if 

+ [u; 3 (r h ) 2 - L0 2 (r h f + if + [3 - LJ 3 (r h ) 2 ] 2 
< 2r 2 h (1 - Arl) . (3.13) 

An example of a typical su(4) EYM black hole was plot- 
ted in [27]]. The solutions have the expected features, 
with the metric functions monotonically interpolating be- 
tween their values on the black hole event horizon and 
at infinity, and the gauge field functions having various 
numbers of zeros outside the event horizon before mono- 
tonically converging to their values at infinity. 

Considering the solution spaces, to produce a two- 
dimensional plot, we now have to fix three parameters. In 
figures [TT] and [TJ] we show examples of the solution space 
when we fix A, and the value of one of the gauge field 
functions at the horizon, varying the values of the other 
two gauge field functions at the horizon. In both figures 
[Til and [T3] we indicate the numbers of zeros of the three 
gauge field functions for the regions where we find black 
hole solutions. Elsewhere in these two figures, we do not 
find black hole solutions. Now that there are three gauge 
field functions, it can be seen that the structure of the 
solution space is quite complicated (and gets ever more 
complicated as |A| decreases). However, for the partic- 
ular values of A and rh in figures [TTJ and [T21 it can be 
seen that there are solutions in which all three gauge field 
functions have no zeros. 




1.6 £0^)1.7 



FIG. 12: Solution space for su(4) EYM black holes with r h = 
1, A = —1 and u>2(rh) = 1-8. Where there are solutions, 
the numbers of zeros of the three gauge field functions are 
indicated for the relevant regions. Elsewhere in the figure we 
find no black hole solutions. As well as the regions indicated, 
we also find small regions where the numbers of zeros of the 
gauge field functions are (711,712,713) = (1, 1,0) and (0, 1, 1). 



Many of the other features of the phase space ob- 
served in the su(2) and su(3) cases are replicated here, 
namely: the fragmentation of the solution space as | A| de- 
creases; as |A| increases, the proportion of the parameter 
space for which the constraint (13.13|) is satisfied and we 
have black hole solutions increases; for sufficiently large 
|A|, we have solutions in which all gauge field functions 
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FIG. 13: Solution space for su(4) EYM black holes with A = 
— 1, th = 1 and uj\{rh) = ^(r^). The shaded region indicates 
those values of the parameters Lu\(rh), L02(rh) for which the 
constraint (|3 . 1 3|) is satisfied, but for which we find no black 
hole solutions. Where we do find black hole solutions, the 
numbers of zeros of the gauge field functions are indicated 
(note that 713 = ri\ in this case). We have also plotted the 
dashed line Auii(rh) 2 = 3uj2(rh)' 2 , on which lie embedded su(2) 
solutions. As well as the regions marked, we also find small 
regions where the numbers of zeros of the gauge field functions 
are ni = 113 = 2, ri2 = and ni — n$ = 0, «2 = 2. 
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FIG. 14: Solution space for su(4) EYM black holes with 
A = —10, rh = 1 and uii(rh) = ^(fh)- The shaded re- 
gion indicates those values of the parameters u\(rh), uJ2(rh) 
for which the constraint (|3.13[) is satisfied, but for which we 
find no black hole solutions. In this case, for all the black hole 
solutions we find, all three gauge field functions have no zeros. 
We have also plotted the dashed line 4aji(Vh) 2 = 3cj2(rh) 2 , on 
which lie embedded su(2) solutions. 

behaviour as is observed at smaller A values in the su(2) 
case. 



have no zeros. Figures [T31 and [JJ] illustrate these fea- 
tures. In both figures [13] and [JJ] we have used the ex- 
ploited the symmetry (|2 . 1 5[) of the field equations and 
set wi(rfc) = ojsirh), although it should be noted, from 
figures [TT1 and [T2l that this does not need to hold (that is, 
although the field equations have the symmetry (|2. 15|) . 
it is not necessary for the solutions to have this symme- 
try). In both figures [TBI and [TTI we have plotted the 
line 4wi(r/j) 2 = 3uJ2(rh) 2 = 4w 3 (r^) 2 , on which lie the 
embedded su(2) solutions (|3.3p . It can be seen in both 
figures that the solution space is symmetric about this 
line, as in the su(3) case. Figure IT"4l shows that it is still 
the case that, for sufficiently large |A|, all the solutions 
we find are such that all three gauge field functions have 
no zeros. 

Comparing figures [51 and !131 we see that the proportion 
of the phase space for which the constraint (|3.13p is sat- 
isfied and we have black hole solutions is rather smaller 
than in the su(3) case. This can be understood from 
the scaling f|3. 5|) required to embed the su(2) solutions 
into su(3) EYM. From (j3~5[) . an su(2) black hole solution 
with cosmological constant A and event horizon radius 
Th is embedded into su(3) as a solution with cosmolog- 
ical constant 4A and event horizon radius r^/2 (since 
A3 = 2 (|3.4p ). and into su(4) as a solution with cosmo- 
logical constant 10A and event horizon radius r^/VlO 
(since A 4 = \/T0 (|3.4p V This scaling means that, for 
larger N , larger A values are needed to find the same 



C. Solitons 

The behaviour of the gauge field functions near the 
origin (|2.34p makes finding numerical soliton solutions 
of the field equations (|2.9l2.12p more complicated than 
finding black hole solutions. We define new variables 
/3i(r), . . . /3jv— i{f) which have the following behaviour 
near the origin: 

P j (r)=b j r j+1 +0(ri+ 2 ), j = 1, . . . , N- 1; (3.14) 

where the bj are the constants in the expansion of oj(r) 
(|2.34p . Therefore the gauge field functions take the form 

N-\ 

Pk(r)v k . (3.15) 

k=l 

For each N, we proceed as follows. Firstly, the normal- 
ized eigenvectors v k of the matrix A4n_i (|2.25p are cal- 
culated. We then have the u>j(r) in terms of the /3fc(r) 
from (|3.15p . The expressions (|3.15p are substituted into 
the field equations (]2 .912 . 1 2[) to give differential equations 
for the /3ft (r). The Yang-Mills equations for the Pk(r) 
will be given explicitly for su(3) below. For the Einstein 
equations, the quantity G (|2.13p becomes 

N-l 
k=l 
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because we have normalized the eigenvectors Vk- The 
quantity pg (|2.10p takes a complicated form in terms of 
the ft(r) (which we do not write here), but is readily 
computed in Maple. Further details of this procedure in 
the su(3) and su(4) cases will be outlined below. 

Many of the features of the solution space for black 
holes are seen also in the soliton solution spaces. In par- 
ticular, the solution space becomes more complicated as 
|A| decreases, eventually reducing to the asymptotically 
flat solution space as A — > 0. As |A| increases, we find 
more solutions and, for sufficiently large |A|, we find so- 
lutions in which all the gauge field functions have no ze- 
ros. In the following subsections, we have focused on the 
structure of the solution spaces for smaller values of |A| 
where there are more features. 




10- 4 io- 3 io- 2 io- 1 10° io 1 io 2 10 3 10" io 5 r io 6 



su(3) solitons 



In su(3) EYM, the matrix M N -i <|2~2"5) with N = 3 
takes the form 



M 2 = 



4 -2 
-2 4 



(3.17) 



It is straightforward to confirm that the eigenvalues of 
M.2 are 2, 6, with corresponding normalized eigenvectors 



v 2 



=-( 1 



(3.18) 



As described above, we therefore write the gauge field 
functions as follows, from (|3.15[) : 



wi(r) 
w 2 (r) 



V2 
V2 



The Yang-Mills equations 
equations for the ft(r): 



V2 \Pi(r)-02{r) J 

then give the following 







2m - 2r J 



2Ar 3 



J (2 + ft) (ft 2 + 4ft + 7ft 2 ) 
2Ar 3 



1 



2m - 2r J 



ft 



fi'2 



- (7ft 2 + 28ft + ft 2 + 24) ft; 



(3.20) 



where, in this case, the expression for pg (|2.10p is not too 
complicated: 



(ft 2 + 4ft + ft 2 ) 



1 

8^ 

+48ft 2 + 48ftft 2 + 12ft 2 ft 2 



(3.21) 



We then numerically integrate the field equations 
(I2.12l3.20p with the initial conditions (|3~T4"]) . The so- 
lution space is described by three parameters: A, ft and 
62. 



FIG. 15: Typical su(3) soliton solution, with A = —0.1, 61 = 
0.35 and b 2 = 0.115. The gauge field function u\{r) has two 
zeros, while u) 2 {r) has no zeros. 
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FIG. 16: Solution space for su(3) solitons with A = —0.1. 
Where we find solutions, the numbers of zeros of the gauge 
field functions are indicated. Elsewhere in the parameter 
space we do not find solutions. For this value of A, there 
is a very small region (near the origin) of solutions in which 
both gauge field functions have no zeros. Although they are 
too small to show on this figure, we also find regions where 
the numbers of zeros of the gauge field functions are n% = 3 
(with n 2 G (0, 1, 2, 3)) or n 2 = 3 (with m E (0, 1, 2, 3)). 



A typical soliton solution is shown in figure 1151 In 
figure 1151 we have plotted the auxiliary functions ft (r) 
and ft(r) as well as the physical field quantities ra(r), 
S(r), u>\{r) and u>2{r). It can be seen that all variables 
have the expected behaviour, both near the origin and 
at infinity. At infinity, the ft(r) functions converge to 
constant values, which can be arbitrary (since the values 
of the gauge field functions u>j{r) at infinity (|2.20[) are 
arbitrary) . 

The solution spaces for two particular values of A are 
shown in figures [TB] and 1171 The origin ft = = 62 
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FIG. 17: Solution space for su(3) solitons with A = —1. 
Where we find solutions, the numbers of zeros of the gauge 
field functions are indicated. Elsewhere in the parameter 
space we do not find solutions. We now have a much larger 
region of solutions in which both gauge field functions have 
no zeros. 



corresponds to pure adS space. In both figures [16] and 
I17[ we see that the solution space is symmetric about the 
axis 62 = 0. This is due to the symmetry (|2.15p of the 
field equations, since the mapping b 2 — * — &2 effectively 
swaps wi(r) and uJiiy) from ()3.19|) . In both figures we see 
a region of solutions in which both gauge field functions 
have no zeros, but the size of this region increases for 
larger |A|. 



2. su(4) solitons 



In su(4) EYM, the matrix M N -i (p?^) reads 

A1, = j -\ 1 2 s - v .12 1 . (3.22) 



6 -VT2 
-7l2 8 -yT2 
-yT2 6 



and has eigenvalues 2, 6 and 12 (|2.28p and normalized 
eigenvectors 



L(f 



1 



1 



1 



1 



(3.23) 



The gauge field functions therefore take the form (|3 . 1 5[) : 

(3.24) 











^3 / 





FIG. 18: Typical su(4) soliton solution, with A = — 10, 61 = 
— 1.2, 62 = —0.1 and bz = —0.01. We plot just the auxiliary 
functions Pk(j') and the gauge field functions i^j(r). In this 
example, all three gauge field functions have no zeros. Note 
that Pz{r) is not identically zero, it monotonically increases 
to 6.41 x 10~ 2 asr^oo. 



The expressions for the Yang-Mills equations (|2.9p and pg 
(|2.10p are now quite lengthy and so we do not reproduce 
them here. We now have a four-dimensional parameter 
space: A, b 1: b 2 and 63. 

A typical su(4) soliton solution is shown in figurc[T51 In 
figure [18] we have plotted just the gauge field functions 
Wj(r) and the auxiliary functions /3fe(^)i as the metric 
functions have similar behaviour to that seen in, for ex- 
ample, figure [T5l In figure [T51 all three gauge field func- 
tions have no zeros. For this large value of |A|, we find 
soliton solutions for a wide range of values of b\ , all with 
the three gauge field functions having no zeros. However, 
we also find that 62 and 63 have rather smaller ranges over 
which we find solutions. This can be seen in the next two 
figures. 

In figures \T§\ and [301 we show the solution spaces of 
solitons for A = —1, for b% = and 63 = respectively. 
In each case we have indicated the numbers of zeros of 
the gauge field functions in those regions where we find 
solutions. Elsewhere, no solutions are found. The solu- 
tion space in figure [12] is not symmetric about the axis 
63 = 0, but, in figure [2TJ1 the solution space is symmetric 
about the axis 62 = 0. This is expected from the form of 
the gauge field functions Uj(r) in terms of the auxiliary 
functions /3fc(r) (|3.24p . It will be seen from figures [Inland 
[2D] that, for this value of A, we have many solutions in 
which all three gauge field functions have no zeros. 



IV. CONCLUSIONS 

In this paper we have presented new soliton and hairy 
black hole solutions of su(iV) Einstein- Yang-Mills theory 
with a negative cosmological constant. Our solutions are 
purely magnetic, so that the gauge field functions are in 
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FIG. 19: Solution space for su(4) solitons with A = — 1 and 
62 = 0. Where we find solutions, the numbers of zeros of the 
gauge field functions are indicated for the relevant regions. 
For other values of the parameters b\ and 63, we do not find 
solutions. Note that the solution space is not symmetric about 
the axis 63 = 0. 
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FIG. 20: Solution space for su(4) solitons with A = — 1 and 
63 = 0. Where we find solutions, the numbers of zeros of the 
gauge field functions are indicated for the relevant regions. 
For other values of the parameters b\ and 62, we do not find 
solutions. In contrast to figure 1191 in this case the solution 
space is symmetric about the axis &2 = 0. 



general described by N — 1 functions, giving N — 1 inde- 
pendent degrees of freedom. This gives, in total, N + 1 
parameters (N — 1 from the gauge fields, plus the cosmo- 
logical constant A and the event horizon radius r^, the 
latter being zero for soliton solutions) which characterize 
the solutions. We have developed the formalism for find- 
ing solutions for arbitrary TV, and have discussed in some 
detail the properties of the solution space for N = 3, 4. 
Although the spaces of solutions get progressively more 



complicated as N increases, the key features are those 
also found in the su(2) case, namely 

1. Solutions exist in continuous open subsets of the 
phase space; 

2. As |A| — > 0, the solution space fragments and ap- 
proaches the discrete solution space in asymptoti- 
cally flat space; 

3. For sufficiently large |A|, we have solutions in which 
all N — 1 gauge field functions have no zeros. 

The last item is of particular interest. The existence of 
these solutions, for sufficiently large |A|, can be proved 
analytically [23|]. In the su(2) case, it is known that 
at least some of the solutions in which the gauge field 
function has no zeros, for sufficiently large |A|, are lin- 
early stable, both under spherically symmetric [H, [l9[ 
and non-spherically symmetric [2(| linear perturbations. 
We have seen how su(2) solutions can be embedded into 
su(N) EYM, and the first question is whether those solu- 
tions which are stable as solutions of su(2) EYM remain 
stable when considered as solutions of su(N) EYM. We 
will show in a separate publication that this is indeed 
the case [231 ] . and, furthermore, that there are genuinely 
Sll(N) solutions, in a neighborhood of these embedded 
su(2) solutions, which are also stable under linear, spher- 
ically symmetric perturbations. The analysis is rather 
involved so we do not describe it further here. The ques- 
tion of non-spherically symmetric perturbations remains 
open at this stage. 

Other interesting open questions remain. Firstly, there 
is evidence [28| that solutions to su(oo) exist in adS, at 
least for sufficiently large |A|. The field equations for 
su(oo) are rather different in structure from those for 
Su(N), with the infinite number of ordinary differential 
YM equations being replaced by a partial differential 
equation (28|. Therefore different numerical techniques 
will be required to solve the field equations. However, 
the fact (to be proved in [23j) that there are soliton and 
hairy black hole solutions in su(N) EYM for any N sug- 
gests that non-trivial solutions of su(oo) EYM may in- 
deed exist. This leaves open the interesting possibility of 
giving a black hole infinite amounts of gauge field hair. 
Secondly, we have not examined the question of whether 
there are topological black hole solutions of su(N) EYM, 
generalizing the topological su(2) black holes found in 
[29(, but we anticipate that such solutions exist. All 
k = su(2) EYM topological black holes are known 
to be stable as are at least some of the k = — 1 solu- 
tions [29( , so the stability of any su(N) EYM topological 
black holes would also be of particular interest. Finally, 
there is the question of the implication of our solutions 
for the adS/CFT correspondence [17]. A black hole with 
a particular mass and magnetic charge measured at in- 
finity in adS can now be cither an abclian, magnetically- 
charged, Reissner-Nordstrom-adS black hole or any one 
of a number of su(N) EYM black holes with different 
N. We would expect that, in analogy with the su(2) case 
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[30( , there are su( N) solutions in some super-gravity the- Acknowledgments 
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